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ABSTRACT 


Asteroseismology is a powerful tool that can precisely characterize the mass, radius, 
and other properties of held stars. However, our inability to properly model the near¬ 
surface layers of stars creates a frequency-dependent frequency difference between the 
observed and the modeled frequencies, usually referred to as the “surface term”. This 
surface term can add signihcant errors to the derived stellar properties unless removed 
properly. In this paper we simulate surface terms across a signihcant portion of the 
HR diagram, exploring four different masses (M = 0.8,1.0,1.2, and 1.5 Mq) at hve 
metallicities ([Fe/H] = 0.5,0.0, —0.5, —1.0, and — 1.5) from main sequence to red giants 
for stars with Teg < 6500Rr and explore how well the most common ways of htting and 
rei noving the sur f ace te rm actually perform. We hnd that the two-term model proposed 
by iBall fc GizonI (120141 ) works much better than other models across a large portion of 
the HR diagram, including the red giants, leading us to recommend its use for future 
asteroseismic analyses. 


Subject headings: asteroseismology - stars: oscillations 


1. Introduction 

Asteroseismic analysis of stars is the most robust way of determining stellar properties. As¬ 
teroseismic data allow us to determine the masses and radii of stars to exquisite precision. NASA’s 
Kepler mission has provided asteroseismic data on hundreds of dwarfs and tens of thousands of 
giants. These data have revolutionized our knowledge of the target stars. Achievements include 
determining the properties of exoplanet hosts, and thereby the radii of the planets themselves (e.g., 


Howell et al 


2012 ; 


Borucki et al 


12012 : 


Carter et al. 


2 OI 2 II. and catalogs of the properties of dwarfs 


(jChaplin et al.l l2014l l and giants (jPinsonneault et al.ll2014l l in the Ke pler field. We can also de¬ 
termine fairly prec i se ages for stars, albeit in a model dependent way (jMetcalfe et al.ll2010l . l2012l : 
Mathur et al. 2012 : Silva Aguirre et al. 2015 . etc.). 


Knowledge of stellar properties, such as age, can illuminate research in a number of differ¬ 
ent fields: planet formation and migration, stellar evolution, and the star formation history and 
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chemical evolution of the Milky Way. The sources of age uncertainties for asteroseismology are 
vastly different from other methods. The best determined ages outside of asteroseismology are 
htting color-magnitude diagrams to theoretical isochrones for star clusters. However, this has the 
obvious disadvantage that it cannot be applied to held stars. Fitting isochrones to single held stars 
is plagued with problems, such as the “terminal age b ias”, the tendency for st ars to be system¬ 
atically pulled towards the end of the main sequence (IPont fc Even 120041.1200511. Even the more 


Takeda et al 


20071 ) that avoid the 


sophisticated, Bayesian models ( Jorgensen fc Lindegrenf20051 : 
terminal age bias often come with large uncertainties. This is further complicated by observational 
difficulties, such as reddening and unknown distances. Asteroseismology, on the other hand, is 
unaffected by reddening or distance. 


For most stars, we only have two pieces of asteroseismic data: the large separation, Aiy, and 
the frequency of maximum power, Additionally, for many of the observed stars, precise aster¬ 

oseismic data, i.e., frequencies of individual modes, have also been observed. Detailed asteroseismic 
analyses can be done for these stars. The usual way to determine the properties of these stars is 
to construct models that reproduce not just the “classical” properties of the star (i.e.. Teg, metal- 
licity, and position on the color-magnitude diagram), but also the pulsation frequencies. Since the 
frequencies are usually extremely precise, the stellar properties can also be determined precisely. 


There is, however, one major issue that plagues asteroseismic analyses, and that is the so- 
called “surface term.” Despite improvements in stellar models, we are still unable to model the 
near-surface layers of a star properly. The major source of uncertainty is how we model convection. 
Convection is an inherently three-dimensional, dynamical phenomenon, and we approximate that 
to a static, one-dimensional e quation either through mix ing length formalism or through other 
approximations, such as that of ICanuto fc Mazzitellil (119911 1. The mixing length approximation also 
ignores the dynamical effects of convection, such as turbulent pressure support and kinetic energy 
flux. Most of the effect of the lack of proper modeling lies very close to the surface of the model, 
notably in the so-called super-adiabatic layer, where convection is inefficient. The consequence of 
this is a frequency-dependent frequency difference between frequencies of the model and that of 
the star, which was first observed in the cas e of the Sun. This is usually called the “surface term” 
(jChristensen-Dalsgaard fc Berthomiedll99lh. The surface term also depends on the atmospheric 
physics used to model the stars. IChristensen-Dalsgaard fc Berthomied (|l99ll ) showed that if one 
corrected for the mode-inertia, then the frequency differences due to the imperfect modeling of the 
near-surface layers is a function of frequency alone. 


Helioseismic analysis of the Sun proceeds through the removal of the surface term using a subset 
of the solar frequencies. This is possible because many thousands of solar oscillation frequencies 
with degrees ranging from I = 0 (the radial modes) to about I = 300 have been precisely determined. 
Filtering out the surface term in the way it is done for the Sun is not possible for other stars; this 
is a consequence of the fact that only low-degree mode-frequencies (/ = 0,1, 2 and sometimes I = 3) 
are observed for these stars. This does not give us enough information to filter out the surface term 
and determine other properties of the star at the same time. In the solar case, we also have the 































- 3 - 


advantage that the mass, radius, and age are known independently, and thus, the seismic analysis 
can focus on determining the details of solar structure. For the vast majority of other stars, the 
mass and radius of a star have to be determined from the same data set that is used to determine 
age and other properties of the star, and thus, filtering out the surface term without putting any 
constraints on what we expect becomes impossible. The usual way out of this bind is to assume that 
the surface term present between the frequency differences of a star and its model can be assumed 
to be a scaled version of the surface term present in the frequency differences between the Sun and 
solar models. The problem with this, however, is that no one has tested if such scaling works for 
all models at different stages of evolution and with non-solar metallicities. Improper removal of the 
surface term can affect the inference derived from th e surface term. They can affect determined 
ages and internal structure (ISilva Aguirre et ahl 1201311 . estimates of the edges of convection zones 


le.g.. lBasu et al.lll994l : iMonteiro et ahl 


1994; 


Basu 


the helium ionization zones (jBasu et al.ll2004l : 


19971: 


Mazumdar et al 


Mazumdar et al 


2 OI 2 I') and the position of 


2 OI 2 I : IVerma et al.1 l2014l l . 


In this paper, we examine whether the usual ways of correcting for the surface term can be 
applied to all stars. For this, we simulate different forms of the surface term between pairs of 
models and determine how well we can remove those differences. The rest of the paper is organized 
as follows. We give examples of the solar surface term in Section [2l The usual ways of determining 
the surface term are described in Section [3j In Section 01 we describe the different ways in which 
we have simulated the surface term. Our results are described in Section [5l 


2. The Solar Surface Term 


As mentioned earlier, surface terms are usually corrected by scaling the solar surface term, and 
thus, one cannot begin a discussion of the surface term without studying it in the Sun. In Fig. H^a) 
we show the raw frequency differe nces between the Sun and the solar model known as Model S 
(jChristensen-Dalsgaard et al.lll99(ll ). and in Fig. lUb), we show the same after the differences have 
been scaled by the quantity: 


Qnl — 


Ei=o{^y 


( 1 ) 


where Eni is the inertia of a mode of order n and degree I, and Ei=q(u) is that of an / = 0 mode of 
at the same frequency as the model with order n and degree 1. This factor effectively corrects for 
the fact that a given perturbation changes the mode frequencies of low-inertia modes differently 
from that of high-inertia modes. Qni = 1 for alH = 0 modes by definition. Note that the scaled 
frequency differences are predominantly a function of frequency. The spread of the points is caused 
by differences in the interior, while the frequency dependence of the difference is entirely from 
near-surface errors. This frequency-dependent frequency difference is the surface term. 


For all low-degree p-modes that are observed in other stars, Qni — 1- Mixed modes, on the 
other hand, can have much higher values of Qni- Since there can be no / = 0 mixed modes, it 
is often useful to define the surface term. Judging by Fig. [H the surface term is quite simple: 0 
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at very low freq uencies and then a st eep change which can be modeled as a simple function, such 
as a power law ([Kieldsen et alJ 120081 1. Unfortunately, however, the surface term is not universal 
and depends on the physics used. This is demonstrated in Fig. [2l where we show the frequency 
differences between the Sun and a number of published solar models. Not only do the surface terms 
look very different, the magnitude of the frequency difference does not tell us anything about how 
different the structure of the model is from the Sun (see Fig. [3]). As is clearly seen. Model S clearly 
has a larger surface term than model BSB, yet in terms of structure of the interior. Model S is 
clearly better. Thus the surface term can mislead us about how well a model matches the star. 

Since multiple physics inputs in the models shown in Figs. [2] and [3] are different, it is difficult 
to pinpoint the exact cause of the differences in the surface term. However, one can show that the 

surface term changes if the model atmosphere is changed, and/or ne ar-surface opacities are changed. _ 

and/ or the formulation of how convection is modeled is changed (see lChristensen-Dalsgaard &: Thompson 
19971 1. In Fig. |3]we show the frequency differences between two solar models, one constructed with 
the Eddington T-r rel ation in the atmosphere and the other with the Krishna Swamy relationship 
(jKrishna Swamvlll966l ). Note that the frequency differences show the typical surface term pattern. 

The structural differences between the two models are concentrated at the surface. 


Thus, to reiterate, the surface term, even that between the Sun and its models, is not unique. 
Thus, even if we assume that the surface term of non-solar models can be corrected using the solar 
surface term, we first need to ensure that the models are constructed with identical physics. 


3. Usual ways of correcting for the Surface Term 


The most common w ay to correct f o r the surface term for stellar models is that proposed 
by iKieldsen et al.l (j2008l l . iKieldsen et al.l (j2008l l noted that the steep part of the surface term of 
Model S could be modeled as a power law and proposed that the surface term for other stars be 
corrected as: 


^corv ^mod — ® 


/ ^mod \ 

V l^rei J 


(2) 


where r'mod is the frequency of any given mode of a model, t'corr is the corrected frequency, a is the 
correction at r'mod = t'ref, and b is the slope from fitting the solar frequency differences. To avoid 
extra free parameters, many investigators use = i^max (jMathur et al.ll2012l : lMetcalfe et al.ll2014l . 
etc.). Recall that r'max scales as: 


M/Mq 


^max,© (R/R0)V©eff/7;ff.o) 


( 3 ) 


(e.g. iKjeldsen fc Bedding 1 199f)l : iBedding fc KieldsenI 120031 ). where M, R, and are the stellar 
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mass, radius, and effective temperature. We refer to the model in Equation [2] as the KBCD08 
model. 

Another way to correct for the surface term is to scale the solar surface term directly. This is 
applied as follows: denote the solar surface term as the vq-Si'q relation. Both i/q and Jz/q are then 
scaled to the mass and radius of the stellar model under consideration using the homology scaling: 


^ - 

Auq 

where the large separation of the model can be obtained either from the scaling relation: 


Av I M/Mq 


( 5 ) 


or by calculating the large separation as an average separation from the frequencies themselves. 
The latter is preferred. The resulting {rvQ-rdh'Q) relation is then used to correct the stellar model 
for the surface term. The surface term of the model is then denoted as: 


^corr l^mod — 


( 6 ) 


where 5i^q^sc is fSi'Q at rz^Q = z^mod- We call this the Scaled Solar method, or ScS-1 (the “ 1 ” 
denoting one free parameter). As a variant, we also add a vertical offset b to Eq. [ 6 ] and call that 
case ScS-2. The vertical offse t mimics the effect of the evolution of the phase term in the expression 
for Az/ as a star evolves Isee lWhite et al.ll2012l b 


Since frequency changes caused by a given perturbation depends on the inertia of the mode — 
high inertia modes being perturbed less than low inertia modes — models KBCD08, ScS-1, and 
ScS-2, when applied to non-radial modes, need to be scaled inversely with the Qn,i for that mode. 
This is particularly needed to ensure that mixed modes, which have very large mode inertia and 
often have Qni > 1 , are handled properly. 


Ball fc GizonI (120141 ) have suggested a surface term correction that is completely independent 


of the solar surface term, which makes it very promising. It has not been widely a dopted as yet, 
which makes testing this model more useful for all those who might wish to adopt it. iBall fc Gizon 
(|2014l f have two models of the surface term 5v. They propose: 


5v{v) 


q3 / ^ y 

E/nl \ Z^ac / 


( 7 ) 


where 03 is a coefficient that is determined by fitting, z^ac is the acoustic-cutoff frequency (recall 
that z'max oc z^ac) and hence the same scaling applies), and E^i is the mode inertia. We call this the 
“BG14-1” model. The second model they propose is: 














u 


3 


-1 


( 8 ) 


5v{v) 


Q3 

Enl 


+ — ^ 

Enl V ^ac / 


which we call the “BG14-2” model. In this form there are two parameters that are determined via 
fitting, 03 and a_i. The explicit dependence of expressions on the mode inertia means that that 
the expressions can be applied to non-radial modes without any modification and obviates the need 
to scale using Qn,i- 


4. Simulating the Surface term 
4.1. Calibrated models with different atmospheres 


Given that the surface term is only known in the solar case, testing out the fitting forms requires 
us to simulate a known surf ace term between models and the stars they r epresent. We are guided in 
this by the solar case where IChristensen-Dalsgaard &: Thompson! (jl997l l showed that changing the 
model atmosphere, surface opacities, and/or the mixing length formalism can simulate a surface 
term. Since changes in the model atmosphere produce a surface term between solar models, we use 
this to construct surface terms for other models too. 


The basic set of models was constructed with the Yale stellar evolution c ode, YREG (jPemaraue et al 


ture opacities (llglesias fc Rogerslll99m. and the 


Nuclear reaction rates were from 


for which we used the Formicola et al 


2008l b We us ed the OPAL equation of state ( Rogers fc NavfonovI I 2 OO 2 II . OPAL high tempera- 


Ade. 


Ferguson et al.l (1200.^ ) low temperature opacities. 
berger et al.l (jl998l ). except for the 7)^^0 reaction, 


(j2004l l result. Unlike standard solar models, we did not 


include diffusion or core overshoot. The initial models were constructed with Fddington T-r re¬ 
lationship and the solar calibrated value of a, which is 1.69 in the absence of diffusion. These 
models were constructed for four masses (M = 0.8,1.0,1.2, and 1.5 Mq) and at five metallici- 
ties ([Fe/H] = 0.5,0.0,—0.5,—1.0, and — 1.5), where we assume that [Fe/H] = 0 corresponds to 
ZjX = 0.023 (|Grevesse &: Sauvallll998l L We evolved the models up the giant branch. Models were 
output at many different ages for the purpose of calculating pulsation frequencies. The evolutionary 
tracks from which our models were picked are shown in Fig. [5j For a subset of these models, we 
constructed models with identical mass, radius, age, and luminosity, but with the Krishna Swamy 
(henceforth, KS) T-r relation using the calibration used to construct standard solar models. 


A calibrated solar model is constructed by requiring a 1 Mq model to have the observed 
luminosity and radius at the solar age of 4.57 Gyr. To achieve this, one notes that there are two 
free parameters in the model — the initial helium abundance, Yq, and the mixing length parameter, 
OL. To get a calibrated model, Yq and a are chosen so that the two constraints (on luminosity and 
radius) are satisfied at the solar age. We follow a similar procedure for models of different masses 
and ages. Our “standard” models are the Eddington models described above, and we require the 
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KS models to have the same radius and luminosity at the same age as the Eddington model, and 
this is again done by changing Yq and a. 

The above calibration method, however, works only for a very narrow range of masses, metallic- 
ity, and ages. The main change between Eddington and KS models is the mixing length parameter 
a, which is what is responsible for the surface term. However, there is a very small change in 
To too. For instance, for the Eddington solar model that we constructed, Iq = 0.27572, but for 
the KS model, it is Tq = 0.27566, i.e., a fractional difference of 0.02%; both models have surface 
ZjX = 0.023. Such a small change in Yq does not really change structure, except in stars with 
convective cores; the small difference in the size of the convective core means that the frequency 
differences are no longer, strictly speaking, due to a near-surface change. For low mass stars with 
radiative cores, the difference begins to matter in the upper part of the main sequence where the 
slight difference in the initial metallicity manifests as a difference in the size of the helium core, 
and hence deep-seated difference between the sound speed profiles of the two models. As a result, 
we have only a very restricted set of models for which we could successfully simulate the surface 
term in this manner. These successful models are listed in Table [H 

The frequency differences between these KS and Eddington models essentially represent the 
surface term. We assume the KS models to be the proxy stars and that they have relative frequency 
errors of lO”'^. We used these errors to generate 100 realizations of the frequencies of these models. 
Thus for each model pair we now get 100 realizations of the surface term. Each of these surface- 
term realizations were then fit using the different functions described in Section [3] to see whether 
they can model surface term between these models. 


4.2. Frequency response for a given near-surface differences 

The theory of linear adiabatic stellar pulsations tells us that the frequency differences between 
two models that are not very different from each other can be expressed as: 




I 



6c^ 

nZ 



(9) 


where bi? jc? is the relative difference between the squared sound speed profiles of two models, bpjp 

2 2 

is the relative density difference, and the func tions Kf’^(r) and (r) are known functions and 
can be calculated for a given model (see e.g., lAntia fc Basdll994l ). Implicit in the equation is the 
assumption that the models have the same mass and radius. 


Equation [9] also shows that we can calculate the frequency differences between two models 
by multiplying a given relative sound speed difference and/or a given relative density difference 
with the kernels for a model and integrating the resulting function. This is the path we take in 
order to simulate the surface term over a large part of the HR diagram. To this end, we calculated 
the kernels for all the Eddington models described in Section 14.11 and applied known near-surface 





differences in sound speed to determine the frequency differences resulting from the difference. 

We consider two distinct cases. In the first case, we use an artificial sound speed difference that 
is confined to the near-surface layers (Fig.[ 6 |), assuming no density difference. The profile is basically 
a truncated Gaussian, and hence we call this set the “Gaussian” set. This is essentially a toy model 
since equations of stellar structure imply that a sound-speed difference will be accompanied by a 
density difference. Nonetheless, this gives us a clean signal to try out the fits. 

The second case is more realistic in that we use self-consistent sound-speed and density dif¬ 
ferences. We constructed two solar models, one with the physics described in Section sn and 
the other where the opacities were modified close to the surface as a function of temperature (see 
Fig. [7K). The changes in opacity are not physically motivated; any similar change will work for 
the purpose of simulating the surface term. The frequency difference between the two solar models 
are shown in Fig. EKb), and the relative squared sound speed and density differences are shown in 
Fig. [ZKc). We then find the frequency response to these sound speed and density differences using 
Eq. m We call this the “Opacity set”. 

The frequency differences obtained for both the Gaussian and Opacity sets were then fit by the 
forms described in Section[3]to see which, if any, does well for which models. As before, we assumed 
relative frequency errors of and generated 100 realizations of the surface term for each model 
by adding random errors consistent with the assumed errors. In order to fit the KBGD08 model, we 
first determined the slope b from the frequency differences obtained with Eq. [9] by using the near¬ 
surface structural differences and the kernels of a solar model constructed with the same physics. 
This was done separately for the Gaussian and Opacity sets. The same frequency differences were 
used as the proxy solar surface term to fit the ScS-1 and ScS-2 models. Given the models range 
from main sequence models to red giants, we fit the surface term for different frequency ranges. 
For models with log 5 > 2.7 we fit the simulated surface term in the range t'max i 5Az/. For stars 
with log 5 < 2.7 (i.e., more evolved stars), we fit over the range t'max i 2Az/. 


5. Results 

5.1. Different solar models 


As had been mentioned earlier, the surface term between the Sun and solar models is not 
unique, but rather depends on the model. Thus the first exercise we attempted was to try to fit 
the different models of the surface term to the frequency differences betwee n the S un an d different 
solar models. We used five different solar models: the se ismic model from Antial (1996) which we 
call INVBB, model BSB(GS98) from Bahcall et ah ( 2006), which we r efer to as BSB, M odel S from 


( Ghristensen-Dalsgaard et al.lll996l ). BP04 from iBahcall et al.l ( 20051 ) . and STD from iBasu et al 


( 20001 ) . We only used the I = 0 modes for the hts and assumed a relative error of 1 x 10 ^ on the 
frequencies. We ht the frequency differences in the range r'max i bAv. Note that by definition the 
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models ScS-1 and ScS-2 will fit the data, hence we concentrate only on the other forms. 

Our results are shown in Fig. [8l As can be seen, the fits can be surprisingly poor even for solar 
models. KBCD08 was formulated with the help of Model S, and sure enough, it fits the surface 
term of the model quite well, and the fit can be improved by changing the range of frequencies fit. 
BG14-1 is only able to properly fit the Model S surface term; it fits the surface term for INVBB 
extremely badly. The most promising is BG14-2, which fits the surface term for all the models 
reasonably well, and the fit can be improved by changing the frequency range of the fits. This gives 
us an early indication of the results we can expect for simulated surface term for the non-solar 
models. 


5.2. Calibrated models 


We turn our attention next to the frequency differences between the pairs of models with 
identical properties except for the model atmosphere used to construct them. The fits to the 
frequency differences between the thirteen pairs of models is shown in Fig. [9l The fit is only shown 
for error-free realization of the surface term in each case. For model KBCD08, the power b in Eq. [2] 
was determined from the frequency differences between two solar models that were constructed 
with different model atmospheres (Fig. [4]) and then applied to the other difference s. Since errors 


in so lar oscillation frequencies are very small, a few parts in 10® (see, e.g., Table 1 of lGhaplin et al 


20071 1. we determined b from error-free data. For models ScS-1 and ScS-2, the frequency differences 
shown in Fig. [4] were assumed to represent the solar surface term, and that, in the scaled form, was 
applied to the other differences. Again, the error-free version was used for this. 

In Fig. [TOl we show the goodness of fit in terms of a reduced (Xu) methods. The 

quantity Xu is defined in the conventional manner to quantify the mismatch between the data and 
the model, i.e., 


1 


Xv = 


N — m — I 


N 

E 

2=1 




err 


( 10 ) 


where, 5t'i,surf is the surface term we are trying to fit, is the surface term model (KBGD08, 

ScS-1, ScS-2, BG14-1 or BG14-2) at frequency ui, ai is the uncertainty in Ui, N is the total number 
of points fitted and m the number of free parameters. As usual, the smaller the value of the 
better the fit to the surface term. What we have plotted is the median value of obtained from 
the 100 realizations for each model pair. The vertical “error bars” are the la spread in x^ obtained 
for the 100 realizations. The results are clear. The one-term fitting models, KBCD08, ScS-1, and 
BG14-1, can perform fairly poorly for many of the stellar models. The two two-term fits, on the 
other hand, both perform spectacularly well for these thirteen calibrated models. The extra term 
in the ScS-2 and the BG14-2 fitting methods appear to be needed to improve the fits. 
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5.3. Surface term from known near-snrface differences 

In Fig. dU we show the fits to the surface term for the Gaussian set for the same subset of 
models shown in Fig. [9] and listed in Table [H The one-term fits perform better for the Gaussian 
set than the calibrated set, though again the two-term fits do better, but only marginally so in this 
case. 

As mentioned earlier, we also simulated the surface term for a large number of models, and 
the models shown in Fig. [11] is only a small subset of all the models. The results for the entire set 
of 1321 models is shown in Fig. [12] We have color-coded the points with the median value of xt 
obtained for the 100 realizations of the surface term at given points on the HR diagram. We see 
that most models work reasonably well in this case, the ScS-2 doing marginally better, particularly 
in the red-giant range. 

We show the the results more quantitatively in Fig. [T^ where we have plotted the distribution 
of Xu for fill realizations and all stellar models for each of the five surface-term models. Since it 
is clear the results depend on the evolutionary state of the star, we have divided our sample into 
three bins in log^f. The ScS-2 model appears to work the best in all cases, but the others are not 
too bad either. 

Given that the surface term simulated by determining the frequency response to a sound speed 
perturbation in the form of a truncated Gaussian is not completely realistic, since we assume that 
there is a sound speed difference in the absence of a density difference, we turn our attention to the 
results for the Opacity set. In Fig. [TT] we show the fits to the subset of models listed in Table [T] 
and Fig. [15] shows how the fits behave over the HR diagram. It is clear that, in this more realistic 
surface term, BG14-2 works better than the rest. ScS-2 also works well. Neither ScS-1 nor KBCD08 
are able to fit the frequency differences. BG14-1 works in the main sequence, but not for evolved 
stars. 

In Fig. [16] we show the Xu distribution for the entire sample in the same manner as for the 
Gaussian set. The two two-term models perform the best for main sequence stars. The three one- 
term models struggle. While BG14-1 works in the main sequence, it quickly deteriorates. KBD08 
and ScS-1 so badly everywhere. Both two-term models do well in the main sequence, BG14-2 being 
slightly better. In the sub-giants segment, we start to see BG14-2 pull away from the crowd, while 
the giants undeniably show that BG14-2 is by and far the best surface term model for low logs' 
stars. 


6. Discussion and Conclusion 

The surface term — the frequency differences between a star and its model caused by our 
inability to properly model the near-surface layers — makes asteroseismic analyses difficult. The 
near-surface contribution is usually removed using ad hoc models of what the difference should look 
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like. In this work, we have simulated different forms of the surface term for models that span a 
large part of the HR diagram and tested different models for removing the surface term. 


There are issues with many of the surface-term models even for main sequence stars, and most 
models perform ve ry badly in th e sub- giant and red giant region. However, the recently proposed 
two-term model of iBall fc GizonI (120141 ) works well in all the parts of the HR diagram. The scaled 
solar model with an offset (ScS-2) works reason ably on the ma i n seq uence but not elsewhere. We, 
therefore, conclude that the two-term model of iBall &: GizonI (|2014l l (Eq. [8] of this paper) is the 
model that should be used for future asteroseismic analyses. 


Our results are not affected by the fact that we have chosen a uniform relative frequency errors. 
We have repeated the test on a subset of our models using an error-distribution that has minimum 
errors at t'max and increasing errors away from t'max- Those results are completely consistent with 
our previous results. 

It should be noted surface term corrections can in many cases distort or hide the seismic 
signature from the stellar interior. As a result, stellar models with different properties can fit the 
frequencies of the same star after surface term corrections. The way to way to distinguish these 
models is to use frequency separation ratios. Frequency separation ratios are essentially insensitive 
to the surface term, and hence can be calculated with the uncorrected frequencies of the model. 
There have not been any tests of how models BG14-1 or BG14-2 could distort the signature from 
stellar interiors. To test this we makes use of the model pairs listed in Tabled) Recall that we had 
constructed pairs of models with different atmospheres; the frequency differences between the two 
models is the surface term. We fit this surface term and subtracted it out from the models with the 
KS atmospheres. We then used these ‘corrected’ frequencies to calculate different frequency ratios 
and compared them with the ratios calculated with the original frequencies. In Fig. dZland Fig. dHl 
we show root-mean-square (RMS) differences (in percent) between the ratios calculated with the 
surface-term removed frequencies and those calculated with the original frequencies. We find that 
the BG14-1 and BG14-2 surface-term models do as well, and sometimes even better, than the other 
surface-term models that we have used. Thus, the adoption of this correction poses no more risk 
than adopting any of the others. 

Most of the results, except of course the ones in Fig. [TTland Fig. [181 were obtained using radial 
modes only. The conventional practice has been to define the surface term using radial modes and 
then apply the correction to the non-radial modes after scaling with Qni- The explicit presence 
of the factor E^i in BG14-2 should have, in principle, allowed us to define the surface term using 
both radial and non-radial models. However, for the few cases that we tested, we find that even 
for the BG14-2 model, the best results are still obtained by first defining the surface term using 
radial modes and then applying them to the non-radial ones. This is the case for stars beyond the 
main sequence. For main-sequence stars, it did not matter whether we used just the radial modes 
or all modes. 
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Table 1. Properties of models that could be successfully constructed with two different 

atmospheric models 


Model number 

Mass (Mq) 

Radius (R©) 

Teff (K) 

Age (Gyr) 

[Fe/H] 

1 

0.8 

0.746 

4944 

4.87 

0.0 

2 

0.8 

0.775 

5043 

9.32 

0.0 

3 

0.8 

0.727 

4557 

3.88 

0.5 

4 

0.8 

0.754 

4650 

9.81 

0.5 

5 

0.8 

0.757 

5616 

5.10 

-0.5 

6 

0.8 

0.857 

5800 

10.99 

-0.5 

7 

0.8 

0.765 

6054 

4.57 

-1.0 

8 

0.8 

0.857 

6183 

8.25 

-1.0 

9 

1.0 

0.963 

5754 

2.99 

0.0 

10 

1.0 

1.162 

5852 

7.86 

0.0 

11 

1.0 

0.946 

5325 

3.80 

0.5 

12 

1.0 

1.063 

5422 

8.51 

0.5 

13 

1.2 

1.325 

5812 

3.90 

0.5 
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Fig. 1.— Frequency differences between the Sun (mode set BiSON-13 of iBasu et al.l l2009li and 


Model S of IChristensen-Dalsgaard et al.l (jl996l i. Panel (a) shows the raw frequency differences, and 
Panel (b) the differences scaled by Qni- The points in black are modes of degree I = 0-2. 
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Fig. 2.— Frequency differences between the 1 = 0 modes of the Sun and different solar mod¬ 
els. Model INVBB is the seismic model of Antia ( 1996), model BSB is m odel BSB(GS98) 


of 


Bahcall et al 


Bahcall et al 


(120061 1. Model S is from (IC hristens en-Dalsgaard et al.l Il996l ). BP04 is from 
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Fig. 3.— (a) The relative differences of the squared sound speed between the Sun and two other 
solar models, Model S and model BSB. (b) The relative density differences between the Sun and 
the two models. The results, which were obtained by inverting the frequency differences between 
the Sun and these models, are from iBasu et al.l ([20091). 
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Fig. 

T-t 


4.— Frequency differences for I < 3 models for a solar model constru cted with the Eddingto n 
relation and another constructed with the Krishna Swamy T-t relation [Krishna Swamvl (119661 1 . 
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Tefi (K) 

Fig. 5.— The evolutionary tracks from which models were picked. The different line styles refer 
to different masses and different colors to metallicity. We only used models with Tefj < 6500 K (to 
the right of the vertical, black line). 
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r/R* 


Fig. 6.— An artificial sound speed difference that was used to simulate a surface term using Eq. [9j 
The difference is a truncated Gaussian. 
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Fig. 7.— (a) The near-surface opacity modification applied to the solar model, (b) The frequency 
differences between a normal and a modified opacity solar model, (c) The relative squared sound 
speed difference and relative density difference between the two models. These differences were 
used in Eq. [9]to simulate a surface term in different models. 



V (/iHz) V (a^Hz) V (a^Hz) 


Fig. 8.— The surface term between the Sun and different published solar models fit with different 
models of the surface term in the range t'max i 5Az^. In each panel, the points are the data, and 
the lines are the fits. The different colors and symbols indicate different solar models. 
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Fig. 9.— The frequency differences between the calibrated pairs of models fit with the five surface 
term models in the range r'max i 5 Az 2 . The points are the frequency differences, and the lines are 
the best fit. 
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Fig. 10.— The median xt with la spread from 100 realizations of simulated measurement errors 
for the fits shown in Fig. [9l The vertical, gray lines separate the different models. The model 
numbers correspond to the those listed in Tabled) 
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Fig. 11.— The same as Fig. 0 but for the case in which the surface term was simulated by 
determining the frequency response of a sound-speed perturbation in the form of a truncated 
Gaussian (Fig. [^. 
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Fig. 12.— The median Xu of each stellar model from 100 realizations of measurement errors for 
the different surface term models over the HR diagram for the case in which the surface term was 
simulated by determining the frequency difference caused by an artificial sound speed difference in 
the form of a truncated Gaussian (Fig. [6|). The discontinuity in color in the giant branch, most 
easily seen in the BG14-2 panel, is the switch from changing the fitting range from t'max i bAv 
to t'max i 2Az^. BG14-2’s red branch fits also appear to be moderately worse than BG14-l’s. 
However, this is because the extra fitting parameter does not significantly improve the fit, and the 
xt penalizes fitting formulas with more free parameters. 
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Fig. 13.— The distribution of using all 100 realizations of each stellar model, for the different 
surface term models for three ranges of log 5 for the Gaussian set. Note the change in the scale of 
the x-axis for the low log^f panel. 
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Fig. 14.— The same as Fig. El but for the case in which the surface term was simulated using 
sound speed and density differences derived from near-surface changes to opacity shown in Fig. [3 
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Fig. 15.— The same as Fig. m but for the case in which the surface term was simulated by using 
sound speed and density differences derived from near-surface changes to opacity shown in Fig. 0 
Note the the color scale extends to a much higher value of xt than Fig. [12] and is now logarithmic. 
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Fig. 16.— The distribution of using all 100 realizations of each stellar model, for the different 
surface term models for three ranges of logs' for the Opacity set. Note that the x- and y-scales 
have changed from Fig. [13] and particularly that the x-scale is now logarithmic. 
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Fig. 17.— The RMS difference, expressed in percentages, between the separation ratio ro 2 cal¬ 
culated using surface-term removed frequencies and ro 2 calculated using the original frequencies 
of the models listed in Table [U constructed with KS atmospheres. The RMS differences are plot¬ 
ted against model number. The vertical gray lines separate the models. We plot the differences 
for all the different surface-term models that we have tested in this work. For reference, typical 
uncertainties in calculating the ro 2 (and roi and rio) separation ratios are of order a few percent. 



Fig. 18.— The same as Fig. EL but for the frequency ra tios roi and riQ. We have calculated these 


ratios as per the definitions in ISilva Aguirre et al.l (j201ll ). 























